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Abstract. A new dynamical parameter, the /-indicator, is introduced and used in order to distinguish 
between regular and chaotic motion in galactic Hamiltonian systems. Two kinds of galactic potentials 
are used: (i) a global potential, which describes the whole galaxy and (ii) a local potential, which 
is made up of perturbed harmonic oscillators and describes motion near an equilibrium point. The 
new indicator is based on the energies of the separable system along the x,y and z axis. Comparison 
between the outcomes obtained using the new dynamical parameter and other methods, such as the 
maximum Lyapunov Characteristic Exponent (L.C.E), or the S{c) dynamical spectrum, shows that 
the new dynamical indicator gives fast and reliable results concerning the regular or chaotic character 
of the orbits. The new indicator was tested in several Hamiltonian systems of two (2D) degrees and 
three (3D) degrees of freedom. 



1 Introduction 

Knowing whether the orbits of a dynamical system are ordered or chaotic, is fundamental for the 
understanding of the behavior of the Hamiltonian system. In the dissipative case, the distinction is 
easily made, as both types of motion (regular or chaotic) are attracting. In conservative systems, 
however, distinguishing between order and chaos is often a delicate issue (e.g., when the chaotic or 
regular regions are small), especially in Hamiltonian systems with many degrees of freedom, where 
one cannot easily visualize the structure of the phase space. For this reason, it is of great importance 
to have quantities and indicators, in order to determine if an orbit is ordered or chaotic, independent 
of the phase space. 

Over the last years, several methods have been introduced, in order to determine and characterize 
the nature of orbits. One of the classical ways to do to that is by computing the maximal Lyapunov 
Characteristic Exponent (L.C.E) (see Lichtenberg and Lieberman, 1992). The main drawback with 
the L.C.E, is that this method is very time consuming as in most cases requires time periods of order 
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of 10 time units or even more, in order to obtain reliable results. This fact has lead researchers in 

search of new faster methods, in order to be able to distinguish between order and chaos in galactic 
dynamics and celestial mechanics. During the last years a variety of methods have been presented, 
in order to distinguish order from chaos. Among them it is worth mentioning the Fast Lyapunov 
Indicator - FLI (Froeschle and Lega, 2000), the GeneraUzed Alignment Index - GALI (Skokos et al., 
2007), the Smaller Alignment Index - SALI (Skokos, 2001) and the P{f) method used by Karanis 
and Vozikis (2008). Moreover, another interesting and fast method for distinguishing between order 
and chaos in Hamiltonian systems, is the method of the dynamical spectra of stretching numbers. 
Dynamical spectra have been frequently used in Galactic Dynamics, as fast and reliable indicators 
regarding the nature of motion (see Voglis et al., 1999; Caranicolas and Papadopoulos, 2007; Zotos, 
2011a,b). The reader can find more details and theoretical explanations about the shape and behavior 
of the dynamical spectra, in the pioneer works of Froeschle et al. (1993), Voglis and Contopoulos 
(1994) and Contopoulos et al. (1997). 

The aim of the present article is to introduce and use a new dynamical parameter called the /- 
indicator, in order to distinguish order from chaos in a class of Hamiltonian systems with two (2D) 
or three (3D) degrees of freedom, describing global or local motion in galaxies. There are two basic 
reasons for introducing the /-indicator: The first reason is that it is simple and easy to be applied in 
different kinds of dynamical systems. The second reason is that it provides fast and reliable results. 

The article is organized as follows: In Section 2 we present the galactic potentials and give the 
definition of the /-indicator. In Section 3 we present results for the 2D systems. In the same Section 
we present results obtained using the L.C.E and the 5(c) dynamical spectrums in order to check the 
accuracy of the results given by the /-indicator. In Section 4 we study the nature of orbits in the 3D 
systems. A comparison with other dynamical parameters is also made in the same Section. We close 
with Section 5, where a discussion and the conclusions of this research are presented. 

2 The galactic potentials and the /-indicator 

Two galactic potentials will be used in this research. The first is the logarithmic potential 

VG{x,y,z) = -j-ln[:^ -hc^ + ay^ + bz^ + cl\ . (2.1) 

Potential (2.1) describes global motion in a triaxial elliptical galaxy. Here a,b are the flattening pa- 
rameters, Cb is the scale length of the bulge component, while the parameter X << 1 introduces a 
small asymmetry in the system (see Binney and Tremaine, 1987). The parameter stands for the 
consistency of the galactic units. The choice of potential (1) is justified by the fact that triaxilities are 
common in elliptical galaxies (see e.g. Davies and Birkinshaw, 1988; Van Gorkom and Schimino- 
visch, 1997; Bak and Statler, 2000; Hibbard et al, 2001; Young, 2005). 
The second potential is given by the equation 

Vdx^y.z) = y {x^+/+z^) -8 + +ly{xY+y\^+xV)] . {1.1) 

Here CO is the common unperturbed frequency of the oscillations along the x,y and z axis, (3,7 are 
parameters, while £ is the perturbation strength. Potential (2.2) is a 3D perturbed harmonic oscillator 
and describes local motion in the central parts of a galaxy. Such local 2D or 3D potentials, are products 
of expansion of global galactic potentials in a Taylor series near a stable equilibrium point and have 
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been extensively used in order to describe local motion in galaxies (see Deprit and EUpe, 1991; 

Caranicolas, 1994; Elipe and Deprit, 1999; Elipe, 2000; Arribas et al., 2006; Zotos, 2010). 

The outcomes of this research are based on the numerical integration of the equations of motion 

^^_dV{x,y,z) 



y = -- 



dx 

dV{x,y,z) 



dy 



dV{x,y,z) 

Z= ^ , (2.3) 

where V{x,y,z) is one of the VG{x,y,z) or VL{x,y,z) and the dot indicates derivatives with respect to 
time. The corresponding Hamiltonian reads 

H = ^{pl + pl+p^^)+V{x,y,z)=C , (2.4) 

where Px,Py,Pz. are the momenta per unit mass conjugate to x,y and z respectively, while C is the 
numerical value of the Hamiltonian. We use C = E or C = h, for the global or the local potential 
respectively. 

In this paper, we use a system of galactic units, where the unit of length is Ikpc, the unit of mass 
is 2.325 X 10^ Mq and the unit of time is 0.97748 x 10^ yr. The velocity and the angular velocity 
units are 10 km/s and 10 km/s/kpc respectively, while G is equal to unity. The energy unit (per unit 
mass) is 100 (km/s)^. In the above units we use the values: \)o = 15, Cb = 2.5, a = 1.5, b = 1.7, co = 
0.6, P = 0.2, 8 = 0.4, while X and y are treated as parameters. 

All numerical calculations of the present research, are based on the numerical integration of the 
equations of motion (2.3), which was made using a Bulirsh-Stoer integration routine in Fortran 95, 
with double precision in all subroutines. The accuracy of our calculations was checked by the con- 
stancy of the energy integral (2.4), which was conserved up to the twelfth significant figure. We have 
constructed very accurate integration routines, in order to minimize as much as possible the errors 
(noise). For all numerical calculations of the present research, the errors were of order of 10~^^ or 
even smaller and therefore are negligible. 

The /-indicator for a 3D potential is defined as 



where 



E,= ^-pl + V{x,y = 0,z = 0), 
Ey = ^p^y + Vix = 0,y,z = 0), 

Ez=^-pj + V{x = 0,y = 0,z), (2.6) 



are the energies of the separable system (ai = 0), along the x,y and z axis. For the local potential Vl 
we have 
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z 



2 

while for the global potential Vg we take 



\[pI + ^\^)-z^z\ (2.7) 



Ex =\pI + V{xj = 0,z = 0) = \pI + ^In {x^ -iJ + 4) , 
Ey = \py+V{x = 0,y,z = 0) = ^p^. + "^In {a/ + 4) , 

Ez = \pI + V{x = 0,y = 0,z) = \pI + %n {bz^ + 4) . (2.8) 



For a 2D potential the /-indicator is defined as 

Er 

f2 = 4 ■ (2.9) 
Ey 

In order to avoid confusion we would Uke to draw the attention of the reader to the following 
point: the /-indicator is defined by equations (2.5) and (2.9), as if the potentials used were separable. 
But in all cases, Ex, Ey and Ez and consequently the /-indicator, are computed not for the separable 
potentials but for the potentials given by relations (2.1) and (2.2), where the coupling parameters a,b 
and Y are always present and different from zero. This means that Ex, Ey and Ez entering equations 
(2.5) or (2.9), are not any more integrals of motion in all presented calculations. On this basis, there 
is a coupling of the energies of the test particle along the x,y and z direction. This coupling, although 
hidden, is entering in the formulae given by equations (2.5) and (2.9). It is this coupling that defines 
the regular or chaotic character of orbit. 

A second interesting point is that the above method can be appUed to a wide range of potentials, 
such as potential made up of harmonics oscillators or global mass models and logarithmic galactic 
potentials (see Binney and Tremaine, 1987, pp 42). The author would like to make clear that there 
might be cases of more complex potentials such as those describing spiral arms or binary systems, 
where the method of the /-indicator is not easy to be applied. 

Before closing this Section, we would like to recall the definition of the 5(c) spectrum. The pa- 
rameter c is defined as 

~ Pxi ,r, 1 n\ 

Ci = , (2.10) 

Pyi 

where {xi,pxi,Pyi) are the successive values {xi,pxi,Pyi) on the Poincare phase plane. The dynamical 
spectrum of the parameter c is its distribution function 

where AA'^ (c) are the number of the parameters c in the interval (c, c + Ac) , after A'^ iterations. The 5(c) 
spectrum can be used both in 2D and 3D potentials. Note that, the coupling of the third component z, 
is hidden (see also Caranicolas and Zotos, 2010), but in any case it affects the values of x,px and py 
entering equation (2.10). 
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3 Results for the 2D systems 

Let us now proceed to study the character of motion in the corresponding 2D potentials. First we shall 
present results for the global 2D potential 

VG{x,y) = ^ln[x^-'hc^+ay^ + cl] . (3.1) 
The corresponding Hamiltonian to potential (3.1) writes 

HG2 = \{pl + P^y)+VG{x,y)=E2 , (3.2) 

where E2 is the numerical value of the Hamiltonian. 

A simple qualitative way of revealing the dynamical structure of a Hamiltonian system, is by plot- 
ting the successive intersections of the orbits with a Poincare surface of section (P.S.S). This method, 
has been extensively applied to Hamiltonian systems with two degrees of freedom, as in these systems 
the P.S.S is a two dimensional plane. In dynamical systems with three degrees of freedom, however, 
the P.S.S is four dimensional and the behavior of the orbits cannot be easily visualized. One way 
to overcome this drawback, is to project the P.S.S to space with lower dimensions (Caranicolas and 
Zotos, 2009). However, even these projections are often very complicated and difficult to interpret. 

Figure la shows the {x,px),{y = 0,py> 0) phase plane, when X = 0.01 and E2 = 487. As one can 
observe, a large part of the phase plane is covered by chaotic orbits. The regular regions are confined 
around the two stable periodic points on the x axis. These points represent two identical stable periodic 
orbits around the origin traversed in opposite directions, characteristics of the 1 : 1 resonance. There are 
also regular regions in the outer parts of the phase plane. The corresponding orbits are quasi periodic 
orbits characteristics of the 2:3 resonance and box orbits. There are also some small islands produced 
by secondary resonances. These secondary resonances represent complicated orbits of the dynamical 
system (see Fig. 3a). Figure lb is similar to Fig. la but when X = 0.02 and E2 = 494. Here the chaotic 
region is larger. The regular regions around the periodic points on the x axis have decreased, while 
there are not any box orbits. Secondary resonances are also present. 
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The local 2D potential is 

VL(x,y) = y(x2+/)_e[p(x^+/)+2Yxy] , (3.3) 
while the corresponding Hamiltonian is 

HL2 = l{pl+pj)+VL{x,y)=h2 , (3.4) 

where h2 is the numerical value of the Hamiltonian. The results obtained from the study of the 2D 
local system, will be used in order to help us to understand and reveal the structure of the more 
complicated phase space of the 3D local system, which will be presented in the following section. 

The corresponding phase planes are shown in Figures 2a-b. The value of h2 is 0.10125 in both 
figures, while y = —0.2 for Fig. 2a and y = —0.8 for Fig. 2b. In Fig. 2a the phase plane is covered by 
regular orbits except for a small chaotic layer near the separatrix, while in Fig. 2b the larger part of the 
phase plane is covered by chaotic orbits. Regular regions are observed near the stable periodic points 
on the X and px axis. The corresponding orbits are quasi periodic characteristics of the 1:1 resonance. 




Fig. 2: The {x,px) phase plane for the local potential when h2 = 0.10125. 



In the following we shall study the character of orbits in the 2D potential using three different 
criteria. The first is the L.C.E (see Lichtenberg & Lieberman, 1992), the second is the 5(c) spectrum 
(see Caranicolas and Papadopoulos, 2007; Caranicolas and Zotos, 2010) and the third is the new 
/-indicator. Figure 3a-d shows results for an orbit in the 2D global potential (3.1). This orbit is a 
characteristic example of the 3:4 resonance. The initial conditions are: xq = 7.47, /j^o = 0, while for 
all orbits yo = and the value of pyQ is always found using the Hamiltonian (3.2). The values of all 
other parameters are as in Fig. la. In Fig. 3a we observe a quasi periodic. One can see in Fig. 3b that 
the L.C.E vanishes as expected, the S{c) dynamical spectrum produces four U type spectra showing 
that the orbit produces four islands on the phase plane. Looking at the plot of the evolution of the 
/2-indicator with time shown in Fig. 3d, we observe a quasi periodic profile with almost symmetric 
peaks. This suggests, that the corresponding orbit is regular. The integration time of orbit shown in 
Fig. 3a, is 100 time units, while for the 5(c) spectrum is 2 x 10'* time units. 
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Quite different are the results shown in Figure 4a-d. The orbit shown in Fig. 4a is the orbit produc- 
ing the separatrix shown in the outer part of Fig. la. The initial conditions are: xq = —0.5l,pM = 23. 
The values of all other parameters are as in Fig. la. The L.C.E in Fig. 4b has a value of about 0.017, 
while the S{c) spectrum shown in Fig. 4c is asymmetric with a number of small and large peaks, 
which is also an indication of chaotic motion. Finally, the /2-indicator in Fig. 4d appears to be highly 
asymmetric showing abrupt changes. Furthermore, one observes large deviations between the maxima 
and also large deviations between the minima in the [/2 — t] plot, shown in Fig. 4d. These character- 
istics suggest chaotic motion. The above results, strongly suggest that the /2-indicator gives reliable 
results for borderhne chaotic orbits, in the global potential. The integration time of the orbit and the 
5(c) spectrum shown in Figs. 4a and 4c is 200 time units and 2x10^ respectively. 

The next step is to present results for the 2D local potential (3.3). Figure 5a-d is similar to Fig. 
3a-d but for an orbit in potential (3.3), with initial conditions: xq = 0,pxo = 0.3. The value of pyQ is 
always found using the Hamiltonian (3.4). The values of all other parameters are as in Fig. 2a. All 
three indicators show that the orbit is regular. The integration time of the orbit and the 5(c) spectrum 
shown in Figs. 5a and 5c are same as in Figs. 3a and 3d respectively. On the other hand, the results 
given by three indicators in Figure 6a-d, suggest that the corresponding orbit is chaotic. The initial 
conditions of the above orbit are: xq = 0.47, p^o = 0, while the values of all other parameters are as 
in Fig. 2b. As this orbit is a chaotic orbit, we conclude, that the /2-indicator gives reliable results for 
chaotic orbits both in the global and local potentials. The integration time of the orbit and the 5(c) 
spectrum is same as in Figs. 4a and 4c. 

The above procedure was done for a large number of orbits (about 1000) in the global and local 
potentials. The outcomes strongly suggest, that the /2-indicator can be safely used in order to detect 
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(a) A quasi periodic orbit in the 2D global potential. (b) The corresponding L.C.E. 




(c) The 5(c) spectrum (d) The corresponding /2-indicator. 



Fig. 4: For a borderline chaotic orbit in the 2D global potential. See text for details. 

the regular or chaotic character of a 2D orbit. In all cases, the results given by the /2-indicator were 
found to coincide with the outcomes given by the L.C.E and the 5(c) spectrum. Here we must note, 
that the method of the /2-indicator is faster than the other two methods, as we need 2000 - 5000 time 
units in order to obtain reliable results, while for the S{c) dynamical spectrum we need an average 
time of 5000 - 20000 time units, in order to obtain reUable results. Furthermore, the S{c) spectrum 
needs the calculation of the phase plane produced by the orbit, while the /2-indicator needs only the 
orbit. On the other hand, it is beyond any doubt, that the /2-indicator is much more faster than the 
L.C.E, which needs time periods of order of 10^ time units, in order to give reliable results. 



4 The character of motion in the 3D models 

In this Section we shall investigate the regular or chaotic nature of orbits in the 3D potentials. We shall 

start our investigation from the global potential (2.1). All orbits were started with initial conditions 
{xo,Pxo,zo) ,y() = PzO = 0, where {xo,Pxo), is a point on the phase planes of the 2D system. This point 
lies inside the limiting curve 

^pI + Vg{x)=E2 , (4.1) 

which is the curve containing all the invariant curves of the 2D system. We use E = E2, where E is 
the numerical value of the Hamiltonian for the global 3D potential and the value of pyo, for all orbits, 
is obtained from the Hamiltonian (2.4). The numerical calculations suggest that orbits with initial 
conditions {xo,pxo,zo) ,yo = PzO = 0> such as {xo,Pxo)^ is a point in the chaotic regions of Figure la-b. 
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(a) A regular orbit in the 2D local potential. 




(c) The 5(c) spectrum 
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(a) A quasi periodic orbit in the 2D global potential. 
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(b) The corresponding L.C.E. 
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(d) The corresponding /2-indicator. 
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(b) The corresponding L.C.E. 




(c) The 5(c) spectrum (d) The corresponding /2-indicator. 

Fig. 6: For a chaotic orbit in the 2D local potential. See text for details. 
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(a) A chaotic orbit in the 3D (b) The corresponding L.C.E. 

global potential. 




(a) A chaotic orbit in the 3D (b) The corresponding L.C.E. 

global potential. 




(c) The 5(c) spectrum (d) The corresponding /3-indicator. 



Fig. 8: For a regular resonant orbit in the 3D global potential. See text for details. 
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(a) A chaotic orbit in the 3D global poten- 
tial. 
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(b) The corresponding L.C.E. 
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(c) The S(c) spectrum (d) The corresponding /3-indicator. 

Fig. 9: When zo = 0.6. See text for details.. 



for all permissible values of zo, produce chaotic orbits. On the other hand, orbits with initial conditions 
(•^0)P;cO;^o) ,yo = PzO = 0^ such as {xo,pxo), is a point in the regular regions of Figure la-b, behave in 
quite a different way. Orbits with initial conditions, such as {xo,Pxo), is a point in the regular region 
around the two elliptic periodic points, on the x axis, are regular for values of zo ^ 2.7, while for larger 
values of zo they become chaotic. On the other hand, orbits with initial conditions in all other regular 
regions of Figs, la-b, give regular orbits, for zo < 0.46, while they become chaotic when zo > 0.46. 

Figure 7a-d shows the orbit (a), the L.C.E (b), the S{c) dynamical spectrum (c) and the /3-indicator 
(d), for an orbit in the 3D global potential. The initial conditions are: xq = —l.5,pxo = 15, zo = 0.1, 
while the values of all other parameters are as in Fig. lb. It is important to note that, all three dynamical 
indicators show that the nature of the orbit is chaotic. The integration time for all 3D orbits is 200 
time units, while for the 5(c) spectrum is 2 x 10"^ time units. Figure 8a-d is similar to Fig. 7a-d, for an 
orbit in the same potential with initial conditions: xq = —7.4, = 0, zo = 0. 1 . The values of all other 
parameters are as in Fig. lb. Here we see a quasi periodic orbit, characteristic of the 2:3 resonance. As 
expected, once more, all three indicators coincide, showing that the motion is regular. Figure 9a shows 
an orbit with all initial conditions and parameters same as those used in Fig. 8a, but for zo = 0.6. Note 
that, here the character of the orbit has changed from regular to chaotic. It is obvious, that the chaotic 
nature of this orbit is suggested by all dynamical parameters. 
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(a) A chaotic orbit in the 3D local (b) The corresponding L.C.E. 

potential. 




(c) The 5(c) spectrum (d) The corresponding /3 -indicator 

Fig. 10 



Before closing this Section, we shall present some numerical results for the 3D local potential. 
Fig. lOa-d shows a chaotic orbit (a) and the corresponding indicators (b-d). The initial conditions are: 
xo = 0.48 , pxo = 0, zo = 0. 1 . The values of all other parameters are as in Fig. 2b. Obviously, the motion 
is chaotic. In Figure 1 1 a-d we can see a 3D regular orbit with initial conditions: xq = 0.2, p^o = 0,zo = 
0.1. The values of all other parameters are as in Fig. 2a. In this case, the motion is regular. Finally, 
the orbit presented in Fig. 12a has initial conditions and all other parameters same as those used in 
Fig. 11a, but Zo = 0-6. Here the motion becomes chaotic. Generally speaking, the motion in the 3D 
local potential for orbits with initial conditions {xo,pxo,zo) ,yo = PzQ = 0' such as {xo,Pxo), is a point 
in the regular regions of Figure 2a give regular orbits when zo < 0.26, while for larger values of zo the 
motion becomes chaotic. On the other hand, orbits with initial conditions {xo,pxO,Zo) ,yo = PzO = 0^ 
such as {xo,pm), is a point in the regular regions of Figure 2b are chaotic, when zo ^ 0.17, while for 
Zo < 0.17 they remain regular. 

The above procedure was done for a large number of orbits, both in the global and local 3D 
potentials. For a grid of 100 x 100 equally spaced initial conditions {xo,Pxo) on the {x,px),y = 0,py > 
phase plane and for all permissible values of zo, we need about 4h of CPU time on a Pentium 
IV 2GHz PC. The results strongly suggest, that the /3-indicator can be safely applied in order to 
distinguish between ordered and chaotic motion in 3D Hamiltonian systems. In all cases, the outcomes 
derived by the /3-indicator, were found to coincide with the results given by the L.C.E and the S{c) 
dynamical spectrum. 
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(a) A chaotic orbit in the 3D local (b) The corresponding L.C.E. 

potential. 




(c) The 5(c) spectrum (d) The corresponding /3 -indicator. 

Fig. 11: For a regular orbit of the 3D local potential. See text for details. 



5 Discussion and conclusions 

In the present article, we have studied the regular or chaotic character of orbits in two 3D galactic 
autonomous Hamiltonian systems. The global potential (2.1) and the local potential (2.2). The novelty 
of this research is that in this paper, we have introduced and used a new dynamical parameter, the /- 
indicator, in order to detect the regular or chaotic nature of motion, in different kinds of Hamiltonian 
systems. In general, the nature of an orbit in 2D or 3D dynamical systems, can be identified using 
the /2 or /3 indicator respectively. The character of motion can be visuahzed by the profile of the 
[/ ~ A plot- If motion is ordered, it produces a quasi periodic profile, with almost symmetrical 
peaks. On the contrary, in the case of chaotic motion, the profile produced by the /-indicator, is 
highly asymmetric with large deviations between the maxima and also significant deviations between 
the minima in the [/ — t] plot. The importance of this new dynamical indicator, is that we can obtain 
fast and reliable results regarding the nature of motion, especially in Hamiltonian systems with three 
degrees of freedom, in which the structure of the phase space is complicated and therefore, hard to 
interpret. 

In both galactic systems (2.1) and (2.2) we have started our investigation from the 2D potentials. It 
was observed that the global 2D potential displays considerable chaotic regions together with areas of 
regular motion. Responsible for the chaotic motion is the asymmetric term —Xx^. In the local potential 
the chaotic regions are small when y = —0.2, while for y = —0.8 a large chaotic sea is observed and 
the regular regions are small. Our numerical experiments have shown that the /2-indicator gives fast 
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(a) A chaotic orbit in the 3D local (b) The corresponding L.C.E. 

potential. 




-4 -2 2 4 6 1000 2000 3000 4000 

c t 

(c) The 5'(c) spectrum (d) The corresponding /a -indicator. 



Fig. 12: When zo = 0.6. See text for details. 



and reliable results regarding the nature of motion in both 2D galactic potentials. Comparison with 
the L.C.E and the S{c) dynamical spectrum, justifies the reliabiUty of the /2-indicator. 

The /s -indicator is also a useful tool for distinguishing between order and chaos in 3D dynamical 
systems. Comparison with the L.C.E and the 5(c) spectrum, shows that all three dynamical parameters 
give similar results. This shows that the new indicator can be trusted, in order to distinguish regular 
from chaotic orbits. What is interesting is that the L.C.E needs time periods of order of 10^ time units 
and the 5(c) periods of order 10^ time units, while the /-indicator give rehable results for time periods 
of between 2000-5000 time units. Furthermore, the /-indicator needs comparable time of numerical 
integration as the SALI, FLI and GALI methods, while it is much more faster than the dynamical 
spectra of stretching numbers. 

Before closing, the author would like to make clear that the /-indicator, as defined in Section 
2, was based on the energies along the x,y and z axis ignoring the coupling between the three vari- 
ables x,y and z- But all calculations were made in the 3D potentials (2.1),(2.2) and the 2D potentials 
(3.1),(3.3), where the coupling does exist. It is evident, that the complicated mechanism of this cou- 
phng is responsible for the regular or chaotic nature of orbits, shown very clearly by the /-indicator. 

It is our hope, to be able to apply in the future, the new /-indicator in other kinds of interesting 
potentials and in more complicated dynamical systems, such as those describing binary stellar systems 
(see Caranicolas and Innanen, 2009). With additional theoretical work, we will investigate if it always 
distinguishes correctly between ordered and chaotic motion, in Hamiltonian systems, especially in 
dynamical systems with three degrees of freedom. Therefore, one may conclude, that the /-indicator 
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is a very fast and reliable tool, which can be appUed not only in galactic dynamics, but in non-linear 
dynamics in general. 
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